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We apply the Bose-Hubbard Hamiltonian to a three-well system and show analytically that co- 
herent transport via adiabatic passage (CTAP) of N non-interacting particles across the chain is 
possible. We investigate the effect of detuning the middle well to recover CTAP when on-site inter- 
particle interactions would otherwise disrupt the transport. The case of small interactions is restated 
using first-order perturbation theory to develop criteria for adibaticity that define the regime where 
CTAP is possible. Within this regime we investigate restricting the Hilbert space to the minimum 
necessary basis needed to demonstrate CTAP, which dramatically increases the number of particles 
that can be efficiently considered. Finally, we compare the results of the Bose-Hubbard model to a 
mean-field three-mode Gross-Pitaevskii analysis for the equivalent system. 
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I. INTRODUCTION 

The Bose-Hubbard model provides a useful theoreti- 
cal and experimental platform to study the properties 
of strongly interacting bosonic systems. In particular, 
considerable effort has been applied to the properties 
of interacting bosons in a lattice, especially in the con- 
text of ultracold atom research [TJ [2] . The Bose-Hubbard 
model has also been extensively employed in the investi- 
gation of the properties of ultracold bosonic atoms con- 
fined in a double- well potential [3H9]. Additionally, the 
work on double-well systems has been extended to triple 
well systems (TQHT4] . Despite the relative simplicity of 
the Bose-Hubbard model it provides a rich space of phe- 
nomenology to explore, much of which is inaccessible to 
present analytical techniques due to the growth in size of 
the Hilbert space for a large number of particles. Hence 
there is a need to develop analytically tractable prob- 
lems for Bose-Hubbard systems to gain insight into the 
analytically intractable problems. Here we utilize the 
Bose-Hubbard model to investigate the adiabatic trans- 
port of ultracold bosonic atoms confined in a triple-well 
potential. 

For a single quantum particle in a three-well system 
it is possible to transport the particle between the left 
and right wells such that the probability of finding it at 
any time in the classically accessible state in the mid- 
dle well is negligible. The ideas underpinning such a 
transport mechanism stem from stimulated Raman adi- 
abatic passage (STIRAP) [T5HT9]. In general, STIRAP 
is employed to transfer population between two atomic 
states, via an intermediate state, with the occupation of 
the intermediate state strongly suppressed, and is used in 
quantum optics for coherent internal state transfer [T9l - 
[22] . In analogy with STIRAP, it is possible to adiabat- 
ically transport a quantum particle from the left well to 
the right well such that the occupation of the middle 
well is exponentially suppressed. For such a system the 
three states are typically the ground states of each of the 
three wells and the coupling between the states is the 
tunneling rates between the wells, although it is possi- 



ble to transport between excited states of each well [23] . 
This spatial transport is referred to as coherent tunnel- 
ing via adiabatic passage (CTAP) and has been proposed 
to transport single atoms [24, 25], Cooper pairs [26], spin 
states [27], electrons [23 , 28 - 31, ultracold atoms in atom- 
chip waveguides [32] and photons in three-channel opti- 
cal waveguides [33]. For the case of three-channel op- 
tical waveguides, CTAP for massless particles has been 
demonstrated [34] [35], 

Recently, considerable attention has been paid to the 
CTAP of Bose-Einstein condensates (BECs) in three- 
well systems [36H4Q] . These studies have primarily fo- 
cused on a three-mode model of the system (schemat- 
ically shown in Fig. [TJ, in both the linear and nonlin- 
ear regimes, to determine the robustness of CTAP. Ad- 
ditionally, numerical simulations, using the mean-field 
Gross-Pitaevskii equation, have been performed to model 
CTAP of a three-dimensional BEC containing 2000 inter- 
acting atoms [39]. Here we take a different approach to 
studying adiabatic transport in many-particle systems by 
considering a three-site Bose-Hubbard model with N in- 
teracting bosons. We find that the non-interacting limit 
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FIG. 1. Three- well system where each well is assumed to 
have only one accessible energy level. Wells one and three 
are degenerate, whilst the middle- well ground state is shifted 
relative to the other two by a detuning, A. The tunneling 
rate between wells i and j is . 
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FIG. 2. Hilbert space structure for N particles in three wells. Each layer of the Hilbert space is denoted by N2 and the 
coupling between layers is via K12 (black arrows), which couples sites 1 and 2, and if 23 (red arrows) which couples sites 2 and 
3. 



of this system maps onto CTAP of a single particle along 
a chain of 2N + 1 sites via an alternating coupling pro- 
tocol [4TH43] . This mapping allows us to apply several 
analytical results to the problem, which would otherwise 
be computationally intensive, since the size of the Hilbert 
space basis is of order N 2 . 

This paper is organized as follows: In Section II the 
three-site Bose-Hubbard model is introduced and the 
non-interacting null state is found analytically. The 
breakdown of CTAP in the Bose-Hubbard model is then 
investigated numerically in the presence of interactions. 
In Section III, perturbative methods are used to demon- 
strate how detuning of the middle well can maintain 
CTAP once interactions are introduced into the system. 
Section IV numerically investigates the fidelity of CTAP 
as a function of the interactions and the detuning of the 
middle well. This is compared with an analytic determi- 
nation of the regimes of high fidelity based on the anal- 
ysis presented in Section III. Additionally, a compari- 
son between the results obtained for the Bose-Hubbard 
model and three-mode non-linear mean-field calculations 
is made. Finally, Section V summarizes the key results 
obtained in this work and discusses possible avenues for 
further theoretical and experimental studies. 



define the energy scale of the system by the ground state 
energies of wells 1 and 3, i.e., t\ = 63 = but we include 
a (positive) detuning of well 2 so that 62 = A. 

There are M = (N + l)(JV + 2)/2 ways to distribute TV 
particles amongst three wells and these configurations are 
the Fock basis states |iVi, AT 2 , N3) for the Hilbert space of 
the system. The restriction on the individual well occu- 
pation numbers Ni is Y^=i N% = N. Figure [2] shows the 
Hilbert space structure for the system, with each layer of 
the Hilbert space defined by N2. The Fock basis states 
have a natural ordering based on the individual occupa- 
tions of each well. This ordering gives rise to a triangular 
pattern, with layers defined by the number of particles in 
well 2, and the left to right ordering by the distribution 
of the particles in wells 1 and 3, as depicted in Fig. [2j 
In general, the system can take multiple, looped paths 
through the Hilbert space, giving rise to non-trivial dy- 
namics [45]. However, nonzero detuning on well 2 makes 
higher states difficult to access by the system, reducing 
the contribution from these looped paths. The total adi- 
abatic dynamics can, in the appropriate limit of small 
interactions, be very well described by the lowest two 
layers of the Hilbert space. 



II. THREE WELL BOSE-HUBBARD SYSTEM 

To begin the analysis we consider a three-well sys- 
tem, schematically shown in Fig. [T] subject to the Bose- 
Hubbard Hamiltonian 

H = - K12 (a[a 2 + a\a^j - if 23 (<4^3 + o^a^j 

U 3 3 

+ l^^ni{p>i ~ 1) + 5Z e ^' W 

i=l i=l 

where aj, and hi are the bosonic creation, annihila- 
tion and number operators for site z, Kyi (K23) is the 
tunneling matrix element between sites 1 and 2 (2 and 
3), U is the particle-particle interaction energy and q is 
the ground state energy of site i. We only consider tun- 
neling between nearest neighbors so we set if 13 = 0. We 



A. Non-interacting case 

The Hamiltonian ([!]) cannot, in general, be solved an- 
alytically and must be investigated numerically. How- 
ever, before investigating CTAP of interacting bosons 
in a three-well system, we can gain significant insight 
into the problem by considering the non-interacting case, 
which can, in the adiabatic limit, be solved analytically. 

In the absence of interactions, the system closely re- 
sembles a single particle in a chain of 2N + 1 sites with 
an alternating coupling scheme [15j I4TH43] . It is therefore 
possible to write down the null state of the Hamiltonian 
([!]) for a system of N non-interacting particles in a three- 
well system: 
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FIG. 3. Non- interacting energy eigenvalue spectrum for 
N — 3 using the squared-sinusoidal pulsing protocol Eq. 
(a) With no detuning, A/K = 0, the null state (dashed) is de- 
generate, (b) For small detuning, A/K = 0.2, the states are 
lifted slightly, removing the degeneracy, (c) Moderate detun- 
ing, A/K = 0.6, lifts states enough that some pass through 
the null state and cross at certain times, potentially interrupt- 
ing the adiabatic transport, (d) Large detuning, A/K = 2, 
pushes those crossing states past the null state at the cost of 
bringing the Ei,o- state closer to the null state at t — r/2. 
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The null state is unique for A ^ 0. If the system is 
prepared in the |A))iv eigenstate and evolved adiabat- 
ically, then the system will only occupy states in the 
N2 = layer of the Hilbert space (see Fig. [2J, so there 
is vanishing probability of detecting a particle in well 2. 
The system's Hilbert space trajectory will also explore 
the N2 = 1 layer, since these are necessary to couple 
to the N2 = states, but occupation of the N2 = 1 
states is zero in the adiabatic limit. Hence, under adi- 
abatic evolution, N particles initially in state |iV, 0, 0) 
(if 12 = 0, K23 = 1) can be transported to the state 
|0, 0, N) (K 12 = 1, if 23 = 0), via the states \N - k, 0, k), 
with < k < N. This transformation must be per- 
formed adiabatically to maintain the system in the null 
state, discussed further in Section III. The key issue here 
is that by changing the couplings appropriately we can 
transport the system through the Hilbert space. This is 
equivalent to transporting the particles from well 1 to 
well 3 without intermediate occupation of well 2. 

In the non-interacting limit, the full spectrum of M 
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for integers k and j, where 0<k<N,0<j< [k/2\ and 
Eq = Eq^q = denotes the null state energy. The corre- 
sponding eigenstates \T>kj±)N, hereafter denoted simply 
as \V a V, can be found in closed form for any N. The null 
state (pi exists for any detuning A but, as indicated by 
Eq. (pf7 it is not unique when A = 0, namely when k is 
odd and j = there are eigenenergies that are integer 
multiples of A which vanish when A — » 0. In this case 
there are (N + l)/2 null states for odd TV, and N/2 null 
states for even N. These extra null states have contri- 
butions from basis states with N2 > so if the system 
evolved into one of them then the CTAP transport would 
be disrupted. In addition, the adiabatic coupling (dis- 
cussed in Section III) between the preferred null state 
and these extra null states is undefined, indicating that 
these degenerate states are unavoidable. To prevent this 
we can introduce a small positive detuning to lift these 
degeneracies, as shown in Fig. j3^b). This will isolate the 
desired null state Q and improve the fidelity of CTAP. 
However, larger detunings, shown in Figs. [3|c,d), can in- 
troduce new problems related to the adiabatic couplings 
of other eigenstates. We address this issue in Section III 
B. 

If well 2 is not detuned, i.e. A = 0, then CTAP is still 
possible provided we initialize the system correctly and 
evolve the system adiabatically, despite the extra null 
states [44]. This is because when if 12 = 1, if 23 = or 
vice-versa, the extra null states are in a superposition 
of more than one basis state, whereas the preferred null 
state, Eq. ([2|, is the only eigenstate that is exactly a sin- 
gle basis state. For the case of two particles this is shown 
explicitly below. Hence, we conclude that these extra 
null states do not interfere with CTAP and demonstrate 
this for the case of two particles below [43] . 

The time-evolution of the system is determined by the 
pulsing of the potential barriers K\2 and K 23. If all the 
particles are initialized in well 1, i.e. ^(0) = |iV, 0,0), 
then to induce the required population transfer we pulse 
the left and right barriers in counter-intuitive order. 
The precise form of the pulsing is not as important for 
CTAP as the order, so for convenience we choose squared- 
sinusoidal pulsing 



#12 0) = Ks ' m 



2 'gV K 2 s(t)=Kco^(^), (5) 



V 2r 



where K is the maximum tunneling rate and r is the 
total pulse time. From Eq. ([2|, at t = (t = r), the 
probability of being in the |iV, 0,0) (|0,0, iV)) state is 
1. The probability that the system is in a particular 
basis state \<j>) = |7V l5 N 2l N%) as a function of time is 
= \(N u N 2 ,N 3 \9(i)) | 2 . These probabilities de- 
scribe the trajectory the system takes through the Hilbert 
space over the course of the pulsing. Fig. [4] shows the 
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FIG. 4. Basis state occupation probabilities P^if) for the 
non-interacting null state for ten particles. The evolution is 
for the squared- sinusoidal pulsing, with [7 = and A = 0, in 
the adiabatic limit (r large). The basis states with N2 > 
are negligibly occupied throughout the pulsing. 



population transfer, for ten particles, from the 1 10, 0,0) 
state to the |0,0, 10) state via the null state, Eq. ([2]), in 
the adiabatic limit. 



To demonstrate CTAP in the non-interacting case we 
revisit the case of one particle [24] [25] l28ti30] . The single 
particle basis states are 1 0,1,0), |1,0,0), |0, 0, 1) for which 
the Hamiltonian ([!]) is expressed as the matrix 
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and the null state, Eq. is 
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which has no contribution from |0, 1, 0). This is precisely 
the form of the canonical null state used for CTAP [29] 
and is equivalent to the more familiar dark state in dou- 
bly driven three- level atoms [36] . Similarly, we may con- 
sider the two-particle system, which is spanned by the 
basis states |0,2,0), |1, 1, 0) JO, 1, 1), |2,0,0), |1,0,1) and 
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and the null state, Eq. becomes 

K 2 3 \2, 0, 0) - V2K 12 K 23 \1, 0, 1) + K? 2 |0, 0, 2) 



|2?o> 



(9) 

This state also has no contribution from states with occu- 
pation of well 2 (N 2 > 0). As mentioned above, however, 
if A = there is another null state 
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If the system is in the primary null state, Eq. ([9|, then it 
can be transformed from the state |2, 0, 0) when K\ 2 = 
and if 23 = 1 to the state |0,0,2) when K\ 2 = 1 and 
K 23 = 0, but for both of these cases we see from Eq. (10) 



that the secondary null state \V ) always has a contri- 
bution from the |0,2,0) state. If both particles are ini- 
tialized in well 2 then the system will be in the preferred 
null state \Vq) and will remain in that state if evolved 
adiabatically. 



B. Interacting system 

The value of the Bose-Hubbard model is that it can 
model the interesting physics of real interacting atomic 
systems. The Hamiltonian for the three- well system 
Eq. does not have a general null state for U ^ 0, so to 
find the behavior of the particles we must numerically in- 
tegrate the time-dependent Schrodinger equation for the 
wavefunction ^(t) of the system. Initializing the system 
as *(t = 0) = |vY,0,0), we apply the CTAP pulsing se- 
quence, Eq. (jSp, and determine the probability P| <£)(£) 
that the system is in a particular basis state \<j>) as a 
function of time. If \I/(t) = |0,0,7V), then all the parti- 
cles are said to have been transported across the chain, 
but for efficient CTAP we also require that throughout 
the pulsing protocol the contribution to the wavefunction 
from basis states with N 2 > is minimized. 

In Fig. [5] we plot P l<f)) (t) for r = 500K~\ N = 5 and 
s vam( 



various values of U and A. While five particles is still in 
the few-body regime, rather than the many-body regime, 
it is enough to demonstrate the interesting physics of the 
model while still being numerically tractable with stan- 
dard techniques. Fig. |5|a) shows transport for U = 0, 
A = and the system is confined to the lowest (N 2 = 0) 
layer of the Hilbert space, as expected from the analysis 
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FIG. 5. Hilbert space trajectories P\^ (t) for five-particles, us- 
ing the pulsing sequence denned in Eq. (5), with r = 500/ K. 
(a) U = 0, A = 0, the system evolves as expected, remain- 
ing in the N 2 = subspace. (b) NU/K = 0.03, A = 0, 
strong interactions push the system out of the N 2 = sub- 
space, and reduce the fidelity of the transport even within the 
N 2 = subspace. (c) NU/K = 0.03, A/K = 3, a large de- 
tuning of the middle well restores CTAP and the particles are 
transported across the well with negligible occupation of the 
middle well. In each plot the states are presented along the 
y-axis and are ordered from bottom to top: |0, 5, 0); |0, 4, 1) 
|1,4,0); |0,3,2); |1,3,1); |2,3,0); |0,2,3); |1,2,2); |2,2,1) 
|3,2,0); |0,1,4); |1, 1, 3);|2, 1, 2); |3,1,1); |4,1,0); |0,0,5) 
|1,0,4); |2,0,3); |3,0,2); |4,0,1) and |5,0,0). The dashed 
line marks the desired final state |0, 0, 5). 

presented in the previous section. In Fig. |5|b) the in- 
teractions are increased to NU/K = 0.03, with A = 0, 
and there is a significant departure from the idealized 
null state transport observed in Fig. [5|a). Specifically, 
the transport is no longer confined to the lowest layer 
of the Hilbert space and P|o,o,5)( r ) < 1- However, it is 
possible to recover CTAP via the introduction of detun- 
ing, A. For example, Fig. (5^c), where NU /K = 0.03 and 
A/K = 3, shows full population transfer between wells 



1 and 3, such that P|o,o,5)( r ) ~ 1- However, there is a 
limit to how effective this control is since increasing A 
reduces adiabaticity for constant pulsing time r, as dis- 
cussed below. Increasing the pulsing time r (short of the 
true adiabatic limit r — >> oo) does not restore CTAP when 
A = for an interacting system and can even make it 
worse [45 . The non-adiabaticity in this regime is instead 
due to the CTAP state crossing with other eigenstates, 
similar to Fig. [3jc) for the non-interacting case. 



III. PERTURBATION ANALYSIS 

A. Large detuning limit 

We have seen that making A large relative to NU re- 
stores CTAP for finite pulsing periods, motivating an an- 
alytic perturbation investigation of the Hamiltonian. As 
such, it is beneficial to rewrite the Bose-Hubbard Hamil- 
tonian as H = H/A = (°) + UV, where 



= -K12 (a[a2 + a^i) _ ^23 ^2^3 + a\a2^j 
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K12 = K12 /A, K 2 3 = K 2 $/A and U = U/A is the small 
perturbation parameter if NU/A « 1. In this Section 
the superscript (0) denotes non-interacting quantities. It 
is important to note that although the effect of interac- 
tions changes with time as the distribution of particles in 
the system changes, V is time-independent in the basis 
of Fock states. 

In this limit it is possible to apply standard perturba- 
tion theory to calculate the effect of small interactions 
on the null state of the non-interacting case, Eq. (J2|. To 
first order, the interacting null state energy is 



Eq — E, 



(0) 




"0 



(13) 



and the first-order correction to the null state wavefunc- 
tion is 



m = \^) + ^ { ^PSr \m, (14) 



E (0) 



where a represents the fc, j± of Eq. Q. The eigenstates 
|£>i 0) ) of the non-interacting Hamiltonian corresponding 
to the energies of Eq. @ can be found in closed form, but 
for large N it becomes more efficient to evaluate Eq. (14) 



numerically. From these first order perturbative eigen- 
states it is then possible to calculate the basis state oc- 
cupation probabilities Pi^(t), to investigate the stability 
of CTAP in the presence of weak interactions. In the 
regime of small interactions relative to the detuning of 
well 2 this approach is more efficient than integrating 
the Schrodinger equation. 
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FIG. 6. (a) P|0)(t) numerically calculated for the complete 
Hilbert space (red solid curve) and via perturbation theory 
(dashed black curve) for N = 5, A/K = 1, tK = 500 and 
NU/A = 0.01. In both cases the occupation of states with 
N2 > is too small to be visible, (b) Hilbert space trajectory 
for a two- layer Hilbert space for N = 50, A/K = 1, tK — 500 
and NU/A = 0.01. 

Fig. [6ja) compares the perturbation theory method to 
a complete Hilbert space calculation for five particles. 
For small interactions relative to the detuning, transport 
from the state |5, 0, 0) to |0,0,5) is dominated by the 
7V 2 = states, just as in the non-interacting case. 



B. Adiabaticity 

To determine when it is appropriate to describe trans- 
port of the particles from well 1 to well 3 via CTAP, the 
adiabaticity of the evolution needs to be verified. Increas- 
ing the pulsing time r will in general make the trans- 
port more adiabatic, assuming no eigenstate crossings, 
but otherwise the adiabaticity of the transport within 
the null state \Vq) depends on the coupling to the other 
energy eigenstates \V a ) and their energies E a . To quan- 
tify the adiabaticity of the transport we introduce the 
parameter 



l(2>a|i>0>| 



(15) 



and say that the evolution is adiabatic if A a <C 1 for 
all a 7^ 0. Eq. (15) may be generalized to any pair of 



eigenstates but here we are interested in the coupling 
of the CTAP null state to other nearby eigenstates. It 
is also clear from Eq. (15) that removing degeneracies 



via non-zero detuning A improves the adiabaticity of the 
transport. In particular, we focus on the large detuning 
limit A/K > 1. 

In the absence of interactions, the eigenenergy that 



is 



closest to Eq (for any N) is E 



(o) 
1,0- 



= (A 



y/4(Kf 2 + if| 3 ) + A 2 )/2, from Eq. @. It is then pos- 
sible to determine \V^l_) and calculate the adiabatic 
coupling Ai y o- for the squared-sinusoidal pulsing. Re- 
quiring Ai y o- <C 1 will suppress the coupling between 
this eigenstate and the null state. Since this is the eigen- 
state closest to the null state, A a will most likely be less 
than Ai,o- for all other eigenstates a. From Fig. [3^d) 
we can also see that e[°q_, and therefore Al,o-> is maxi- 



mized at t = r /2 for the pulsing scheme given by Eq. ([5| . 
For A/K > 1 at t = r/2 

* 0- = ^ + ^ + ^£ + O(A-). (16) 



tK 2 



2tA 



16tA 3 



The leading order of Eq. (16) is proportional to A 



so there is a limit on the detuning before the adia- 
batic coupling between the null state and the close state 
|^i°o_) becomes too large and the transport becomes 
non- adiabatic. To first order in A, the requirement 
^1,0- < 1 therefore becomes tK/(2tt\^N) > A/K for 
the non-interacting system. Note that in the true adia- 
batic limit, r — » 00, the adiabaticity parameters A a will 
always vanish. 

In the limit of small interactions relative to the detun- 
ing A ^> NU, we again use perturbation theory to ap- 
proximate the adiabaticity parameter by finding the first- 
order correction to all the eigenstates. Inserting Eqs. ( 13) 
and (|14|) in Eq. (15) we obtain 



A a 



(!)|-n(°)\ 



>(°)|f)( 1 )\ 



E (o) 



E 



(0) 



E (o) 



E 



(0) 



(v^\vT)(eV-e^) + o 



fNU 



where 



|^ 1} ) = E 

07^ 



{vf\V\V^) 



E 



(0) 



(0)\ 



={VW\V\VW), 



(17) 

(18) 
(19) 



which is valid for any pulsing protocol. As with the non- 
interacting case, the adiabaticity of the transport in the 
interacting case can be improved by increasing the puls- 
ing time r, but for finite time, as required in experiment, 
Eq. (17), can be used to show where the CTAP protocol 



becomes non-adiabatic and is used in the next section in 
conjunction with numerical results. 



IV. RESTRICTED HILBERT SPACE 

For the case of small interactions relative to the de- 
tuning, if the system is evolved adiabatically with all 
particles initially in the state |iV, 0,0), the CTAP pro- 
tocol only probes those states which have zero (N2 = 0) 
or singular occupation (N 2 = 1) of well 2, i.e. the bot- 
tom two layers of the Hilbert space structure, shown in 
Fig. [2] If the system is not evolved adiabatically, or the 
interactions are large, then the full Hilbert space con- 
tributes to the dynamical evolution of the system, as 
seen in Fig. |5|b). From this observation it is possible 
to consider the following question. If the Hilbert space 
is restricted to the lower layers (N2 = 0,1,...), does 



the solution in this limit represent a valid approximation 
for the adiabatic transport between states |7V, 0,0) and 
|0,0,7V), via CTAP? Hilbert space restrictions, if valid, 
will facilitate the investigation of systems for large N. 

Rather than look at complete Hilbert space trajecto- 
ries for every different value of A and U to judge the 
validity of a restriction, we can consider two criteria 
to identify successful CTAP. Firstly, to determine if all 
the particles have been transported across the chain we 
calculate the probability that at the end of the pulsing 
P|o,o,iV) ( r ) ~ 1- Secondly, to ensure that occupation of 
well 2 has been suppressed during the protocol we cal- 
culate the maximum occupation of a state with N2 = 1, 
maxP|jy 1) iv 2 =i,iV3}W- The latter can be thought of as a 
lower bound on occupation of well 2, since, if occupa- 
tion of a N2 = 1 state becomes too significant at any 
point, then the system may access higher basis states, 
which reduces the fidelity of the CTAP protocol. These 
two conditions are not mutually exclusive but in tandem 
they can be used to determine the fidelity of the popula- 
tion transfer. Figs. [7^c-f) plots these two quantities for 
a two- and four-layer Hilbert space restriction, compared 
to the full Hilbert space calculation, Fig. (7^a,b) for five 
particles. 

The strongest restriction is to consider only the two 
bottom layers of the Hilbert space structure, states with 
N2 = 0, 1. There are only 2N + 1 basis states so the 
problem scales linearly with iV, rather than with N 2 for 
the full system, and the Hamiltonian can be written in a 
tridiagonal form. The restricted Hamiltonian is still not 
invertible analytically for general U but the Schrodinger 
equation can easily be integrated numerically for a sys- 
tem of several hundreds of particles, since the problem 
is then linear in N (there are 2A^ + 1 basis states in the 
two bottom layers, of the Hilbert space structure). Fur- 
thermore, the null state of this restricted Hamiltonian 
in the non-interacting limit is exactly the null state for 
CTAP of a single particle along a chain of 2N + 1 sites 
[HJ[42]. Fig. [6^b) shows the Hilbert space trajectory for 
N = 50 and NU/A = 0.01 using the two-layer restriction 
(A2 = and N2 = 1). In this figure the only states with 
significant occupation are those with N2 = 0, with es- 
sentially zero probability of the states with N2 = 1 being 
occupied. Since only the states with N2 = 1 couple to the 
rest of the Hilbert space structure (see Fig. [2|, it could 
be expected that in this case a two layer Hilbert space 
is a good representation of the transport that would be 
observed in a full Hilbert space calculation. 

The two-layer restriction, while powerful, is only valid 
in certain regimes. In particular, Fig.[7^c) shows that the 
two-layer restriction allows full population transfer in the 
case of high interaction strength but low detuning, while 
Fig. [TJd) shows that at higher interaction strengths well 
2 is being occupied, regardless of detuning, in contrast 
to the full Hilbert space calculation, Figs. |7^a,b). Away 
from the small interaction limit, these results show that 
the two-layer restriction does not model CTAP well. 

To incorporate larger interactions we can sacrifice ease 
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FIG. 7. (a,c,e) [white=l, black=0] P\ ,o, N )(t = r) for: (a) 
a full Hilbert space calculation; (c) a restricted two-layer 
(N2 = 0, 1) Hilbert space calculation and (e) a restricted 
four-layer (N2 = 0,1,2,3) Hilbert space calculation. (b,d,f) 
[white=0] maxP|jv li jv 2= i ) iV3) (t) for: (b) a full Hilbert space 
calculation [black=0.38], with N2 = 1; (d) a restricted two- 
layer (N2 = 0, 1) Hilbert space calculation [black=0.35], with 
N2 = 1 and (f) a restricted four-layer (N2 = 0, 1, 2, 3) Hilbert 
space calculation [black=0.45], with N2 = 1. In (a,c,e) the 
solid dashed line is given by the energy level criterion, Eq. ( 23 ) 
and solid curves are for Ai,o- = 0.05,0.06,0.07 as calculated 
from Eq. ( |17| ). (g) [white=l, black=0] Meanfield calculation 
of the final population of well three: Ns(t = r). (h) [white=0, 
black=l] Meanfield calculation of the maximum population of 
well two: max N 2 (t)/N. For the above plots tK = 500 and 
for (a-f) N = 5. 

of computation for a larger Hilbert space. A weaker 
restriction is to consider the four bottom layers of the 
Hilbert space tree (A^ = 0,1,2,3). The basis contains 
AN — 2 states, so it still scales linearly with N, but the 
Hamiltonian is not tridiagonal. The benefit is being able 
to investigate couplings to higher states of the Hilbert 
space in the presence of strong interactions. Figs. [7je,f) 
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show the calculations for the four-layer restricted system. 
We immediately see that the inclusion of basis states with 
N2 = 2, 3 is necessary to show that CTAP is unstable in 
the high interaction, low detuning regime, as shown, for 
example, in Fig. [5|b). This concept of expanding the 
Hilbert space from the bottom up, rather than consid- 
ering the whole space, is analogous to the way matrix 
product states minimize entanglement by exploring the 
Hilbert space as required [47]. To determine the point 
at which this instability occurs we need to consider the 
energies of the basis states to find limits in the parameter 
space for when CTAP works. 



A. Limiting the CTAP parameter space 

All the models show that CTAP is not possible in the 
high detuning, high interaction regime. Although we 
have seen that detuning can recover CTAP in the pres- 
ence of interactions, too much detuning will still break 
the adiabaticity of the transport. The solid curves in 
Figs. [7^a,c,e) are contours of constant adiabaticity, cal- 
culated from Eq. (Il7|), for A10- = 0.05,0.06,0.07. For 



a given finite pulsing time r these curves accurately pre- 
dict the limit in the parameter space where adiabatic 
transport is possible. Increasing r will push this bound- 
ary further out and allow stronger interactions relative 
to the detuning of well 2. 

The other region of interest is the low detuning, high 
interaction regime, A < NU. In this regime the two-layer 
Hilbert space calculation, shown in Figs. [7jc,d), signifi- 
cantly deviates from the full and four-layer Hilbert space 
calculations, shown in Figs.[7ja,b) and (e,f), respectively. 
In analogy to the non-interacting case in Fig.[3jc), this is 
where the detuning is low enough for some eigenstates to 
cross the interacting CTAP state; the adiabaticity of the 
transport is broken by these crossing states in the small 
detuning regime, and not necessarily at t = r/2, which 
is typically the time when the adiabaticity parameter is 
maximized. To understand the breakdown of CTAP in 
this region we need to consider the energies of the ba- 
sis states in the presence of interactions. Empirically, 
we find that efficient CTAP only occurs when the maxi- 
mum energy of the N2 = states of the Hilbert space is 
less than the minimum energy of the N2 = 1 states of the 
Hilbert space, as shown in Fig. [8] That is, there is a mini- 
mum detuning which makes N2 = 1 states - which couple 
to higher states - energetically unfavourable, therefore 
resisting any attempt by the system to evolve into higher 
Hilbert space layers. For N particles the highest energy 
basis state of the N2 — layer of the Hilbert space is 
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|3,0,0> |2,1,0> |2,0,1> |1,1,1> |1,0,2> |2,1,0> |0,0,3> 
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(a) 



FIG. 8. Relative energies of basis states in the N2 = and 
iV 2 = 1 subspaces for a system of three particles, (a) With no 
interactions (U — 0) and no detuning (A = 0) these states are 
all degenerate, (b) In the interacting system with no detuning 
(A = 0), the interaction energy (J2i UN i( N i ~ shifts 
each state differently (red shifts). Note that this includes 
|1, 1, 1), which has zero interaction energy, (c) Detuning the 
middle well in the interacting system shifts the N2 = 1 states 
by the same amount (blue arrows) . The dashed line indicates 
that there exists a minimum detuning for which all N2 = 1 
states will have higher energy than all N2 = states. 



the Hilbert space is 



^^,1,^ ) ,1,^) 

U {N - 2\ /iV-4 



2 V 2 J V 2 
U [N\ fN -2 
+ 2 (2 ~ 



(22) 



for an odd and even number of particles respectively. 
Hence, the criterion that the highest energy of the N2 = 
states is less than the lowest energy of the N2 = 1 states 
is given by 



A> 



NU(% 



3 > 

4N y 
3 N 
4N y 



odd AT, 
even N. 



(23) 



and the lowest energy basis state of the N2 = 1 layer of 



These inequalities give the minimum detuning needed to 
efficiently transport N interacting particles across the 
three wells. The region to the right of the dashed line 
plotted in Figs. [7ja,c,e) satisfies Eq. (23) for TV = 5, 
while transport in the region to the left is unstable. 



B. Comparison with Gross-Pitaevskii analysis 

We finally turn our attention to previous studies of 
CTAP in dilute gas BECs [36H10 , which have focused 
on a Gross-Pitaevskii mean-field three-mode model. Fol- 
lowing [39] we consider a one dimensional system, with 
each of the three wells containing a single mode, equiv- 
alent to the Bose-Hubbard system schematically shown 
m Fig.[l] The energy of each mode is the ground state 
energy of each well, U is again the interaction strength 
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(equivalent to g in [39]), A is the detuning of the sec- 
ond mode relative to the left and right modes, and 
describe the wavefunction overlap, and hence tunneling 
rate between modes i and j. In the Gross-Pitaevskii anal- 
ysis the quantities of interest are the populations Ni(t) 
of mode i. Using the same squared-sinusoidal pulsing 
of Eq. ([5]), Figs. [7^g,h) plot the final population of mode 
three, Ns(t = t)/N : and the maximum population of well 
two, maxN2(i)/N, respectively. The regimes of high fi- 
delity CTAP in the three-mode Gross-Pitaevskii analysis 
qualitatively agree with the Bose-Hubbard models. The 
two models are typically applied to very different limits, 
the former being a semi-classical mean-field model for a 
large number of particles, and the latter a linear quan- 
tum model that allows for entanglement. The agreement 
of the two models demonstrates the robustness of the 
CTAP protocol. 



V. CONCLUSIONS 

We have considered adiabatic transport of N particles 
across a three-well Bose-Hubbard system and solved the 
non-interacting case analytically for which we obtained 
the null space and energy eigenspectrum in closed form. 
Detuning well 2 relative to wells 1 and 3 makes the null 
state unique and it is this state in which N particles can 
be adiabatically transported from well 1 to 3, or vice- 
versa, without occupation of well 2. Even without de- 
tuning, CTAP is possible if the system is initialized cor- 
rectly. The non-interacting TV-particle three-well system 
exhibits similar dynamics to adiabatic transport of single 
particles along longer chains of 2N + 1 sites. 

The interacting Bose-Hubbard three-well system can- 
not be solved analytically and there is no general null 
state but the time-dependent Schrodinger equation can 
be integrated numerically for few particles. We have 
shown that for stronger interactions relative to the mid- 
dle well detuning, which disrupt the transport, CTAP 
can be maintained by detuning well 2 enough to isolate 
the CTAP state from other energy eigenstates, but not 



so much as to break the adiabaticity of the transport. 

By considering the weakly interacting case as a pertur- 
bation of the non-interacting case we have defined useful 
limits for when adiabatic transport is possible using a 
squared-sinusoidal pulsing. These limits stem from min- 
imizing the adiabaticity parameters, i.e. the coupling be- 
tween the CTAP state and the closest eigenstate, and 
give good agreement to the limit of detuning and inter- 
action strength for which adiabatic transport is possible 
for a finite pulsing time. 

Having established the detuning-interaction regime for 
which CTAP is possible we have investigated restricting 
the Hilbert space to the minimum number of basis states 
necessary to allow CTAP. For zero or very small interac- 
tions relative to the detuning only states with N2 = or 
1 need to be included in the basis. For interaction ener- 
gies comparable to or greater than the detuning, we must 
include states with N2 > 1, but we can still restrict the 
Hilbert space to the lowest layers to allow CTAP. Apart 
from the adiabaticity criterion Eq. (17), which holds for 



the restricted and unrestricted Hilbert spaces, we have 
obtained a limit on what defines this strongly interact- 
ing regime by comparing the energy of individual basis 
states. In this regime we have seen that adiabatic trans- 
port is not stable. 

Finally, we have compared CTAP in the quantum 
Bose-Hubbard model to CTAP in a three-mode mean- 
field Gross-Pitaevskii approximation. The regime of high 
fidelity of CTAP in both models matches and clearly 
demonstrates the robustness of the CTAP protocol in 
both the quantum and semi-classical mean-field lim- 
its. This mapping between the mean-field approach and 
the Bose-Hubbard model reinforces the use of the semi- 
analytic constraints we have developed to determine the 
validity of CTAP. 
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